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Abstract 

In this paper, we provide exact solutions for the extraction of energy from a rotating black hole 
via both the electromagnetic Poynting flux and matter currents. By appropriate choice of a radially 
independent poloidal function A{0), we find solutions where the dominant outward energy flux is 
along the polar axis, consistent with a jet-like colhmated outflow, but also with a weaker flux of 
energy along the equatorial plane. Unlike all the previously obtained solutions ([l|, [g], [?]), the 
magnetosphere is free of magnetic monopoles everywhere. 
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Jet Formation And Energy Extraction 



I. INTRODUCTION 

In 1977, [l| introduced the force-free, stationary, axisymmetric magnetosphere of the Kerr 
geometry as a possible setting for the extraction of energy from supermassive black holes at 
an astronomical scale. To this day, astrophysicists consider the Blandford-Znajek process as 
the leading mechanism for the observed phenomenon of luminous black holes [for example, 
I4, y, Q. Indeed, it is this mechanism we have focused on as well js, IgUs]. 

In an earlier paper [j], we were sucessful in providing the only known class of exact 
analytic solutions to the equations of force-free electrodynamics in Kerr geometry. Although 
finite everywhere in the magnetosphere, these solutions did not appear to be physically 
realistic. In particular, the currents in the magnetosphere were null vectors fields suggesting 
that the charged particles traveled at the speed of light. Additionally, the current vector 
field was inward pointing resulting in an influx of electromagnetic energy. This seemed to 
suggest that the solutions were not physically interesting. 

In this paper, we provide a clear, physically realistic interpretation of the current density 
vector field. In particular, we decompose the infalling null vector field to currents each of 
which have future pointing timelike velocities as these are candidates for physically realistic 
currents. In this decomposition, one of the currents is outgoing. This current provides a 
concrete mechanism for jet formation in black holes. However, the electromagnetic flux 
continues to be inward pointing. To remedy this, we extend the results via a symmetry 
transformation. Briefly, it turned out that the net current vector is proportional to the 
infalling principle null geodesic of the Kerr geometry. Since the Kerr geometry is of Petrov 
type D 2(], it is only natural to ask whether the equations of electrodynamics would allow 
the existence of solutions where the current vector was proportional to the second principle 
(outgoing) null geodesic of the Kerr geometry. As we shall show in this paper, this is indeed 
the case. And the existence of this dual class of solutions is a general property of the 
equations, not necessarily dependent on our particular solution. The dual solution to our 
particular solution (hereby referred to as the Q_ solution) does allow for the extraction of 
energy via the electromagnetic Poynting flux. 

All generalities are restricted by picking a concrete example and carrying out the energy 
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extraction rates from matter currents and the electromagnetic Poynting flux. Here we flnd 
that the matter current and the electromagnetic Poynting flux naturally describe a polar 
jet for a speciflc choice of an arbitrary poloidal function A. This solution also describes a 
secondary local maxima in the energy extraction rate peaking near the equatorial plane, 
which might correspond to outflowing Poynting flux that could drive a disk wind. Even 
though our past f2+ solution ^ generalized the Blandford-Znajek split monopole solution, 
the non-existence of a magnetic monopole for our new solution is shown here by direct 
computation [this was not the case in the original approximate solution presented byll|. 



II. THE SOLUTION 

In j?] we derived the following class of exact solutions for the force- free magnetosphere 
of the Kerr black hole. Here, the components of the electromagnetic flelds in the Boyer- 
Lindquist coordinate system are given by 
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Please see Appendices |X] and [B] for the deflnitions of the quantities listed above. The general 
subscripts 1, 2, 3 in the Maxwell tensor corresponds to r, 6, ip respectively in our case. Here 
A is an arbitrary function of 6. 

It is only natural to expect the net current vector fleld to follow a geodesic under force-free 
conditions. This is indeed the case: 
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Here n is the infalling principle null geodesic of the Kerr geometry. Explicitly, 

n = — dt - dr + ^ d^. (8) 

A simple calculation will reveal that the solution presented above satisfies Maxwell's equa- 
tions (eq. flBl|) ) and the force- free condition, eq.f lB12|) . The above solution is well defined 
everywhere in the magnetosphere. In particular, since our solution satisfies the Znajek 
regularity condition (eq.( 1B15]) ). the fields are well defined at the event horizon, as we will 
explicitly verify by going into the Kerr-Schild coordinate system (see subsection III A\} . Also, 
the apparent singularity at the poles is removed by the transformation A — )■ sin^ 9 A. 

However, the above solution, as it stands, lacks any meaningful physical interpretation. 
Physically realistic charges cannot fiow along null geodesies. This problem will be remedied 
in the remainder of this paper by decomposing the null current vector into timelike vector 
fields that are possible worldlines of charged particles in the magnetosphere. In Q], it was 
the deduction of a viable expression for Q that immediately gave us the expressions for all 
the fields and currents. Therefore, we shall refer to the complete solution listed above as 
the f2_ solution. 



A. The Solution In The Kerr-Schild Coordinate System 



Transforming the Maxwell tensor F^j, into the Kerr-Schild coordinate system (see section 
lAl). we see that 



and 
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Thus we see that the fields and currents are well defined on the event horizon r = as 
well. This is necessarily so since we had insisted on the Znajek regularity condition given 
by eq.( ]B15l) in the derivation of our solution 7|. 
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B. The Electromagnetic Poynting Flux 

[l| computed the expression for the energy extracted from a rotating Kerr black hole via 
the electromagnetic Poynting vector for the force-free, stationary, axis-symmetric magneto- 
sphere. In the 3+1 notation, the rate of electromagnetic energy extraction becomes 

'^^'''^ ^ H^VtB'^rdA. (12) 



dt 

For the case of our n_ solution, the above expressions gives that 

dSEM Stt r cos2 e 



de<0 . (13) 
dt Jo sin^ 9 - ^ ^ 

I.e., if stationary electromagnetic fields can indeed transfer energy to and from a black hole, 
the solution in particular does not allow for energy extraction. Instead, the black hole 
behaves as an energy sink. We will come back to this point when we consider a dual solution. 

III. A TIMELIKE DECOMPOSITION OF THE NULL CURRENT 
A. Region I 

Region I of the Kerr black hole is defined by the condition r > r+, and the region of 
spacetime given by r_ < r < r_|_ will be referred to as region II. Define vector fields Vj and 
Wi in region I by 

(r^ + g^) dt + ad^ 

and 

^ ^ {3l-l)[{r^ + a')d,+ ad^] + Adr 
^ ^riFA^?>l{?>l - 2) 

Here, / = l{6) is the energy coUimation factor such that l{9) > 2/3 everywhere. Then, 

9{Vi,Vi) = -l and Mm g{Vj, dt) = -1 , 

r— >oo 

and 

g{Wi,Wi) = -l 
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and 



lim g{Wi,dt) 
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That is, Vi and Wj are future pointing timelike in region I, and are candidate proper velocities 
of charged particles. Therefore, in region I we can write the current vector as the flow of 
two oppositely charged timelike currents 
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Naturally, In, describes an outgoing current. For completeness note that in Kerr-Schild 
coordinates 



Vr 



and 
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It is important to note, that in our decomposition, even in the Kerr-Schild coordinate system, 
Vi and Wi are not defined well defined at r+. There is a deeper reason why this happens 
regardless of the decomposition if we require an outgoing current in region I. Consider an 
arbitrary decomposition of the net current vector into terms proportional to future pointing 
timelike vectors in region I. At least one of these vectors, say T must be of the form: 

T = a: + 9^ , 

where is an arbitrary function of spacetime coordinates and X has no other radial com- 
ponents so that the net radial component is positive. The only future pointing causal vector 
at r+ such that its radial component is greater than or equal to zero is of the form 

V+^(f [{rl + a'')di + ad^] 
6 



where d G i? is a constant [see|9|. Then we have that 
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for some non zero constant d (since a unit timehke vector field T cannot abruptly become 
the vector). I.e., T is null at the event horizon (since V+ is null), and hence cannot be 
normalized. This is the reason why V/ and Wi are not defined well defined at r_|_. 



B. Region II 

Here, since we are including our decomposition must be valid at the horizon as 
well. Consequently, all expressions will be given in the Kerr-Schild coordinate system. 
Additionally, the radial component of all the vectors in this decomposition must be inward 
pointing to agree with the causality conditions of the interior geometry (unlike in region I 
where we wanted an outfiowing current). As we shall see, by construction, the decomposition 
here will be well defined at the horizon. In region II define vector fields Vn and Wn by 

[{r^ + a')dt + ad^ + {/\-2)d,] 
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Then, 



and 
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L l^" + ci") 9t + ad^ + {A-l) dr 
Wn = — 



Vp'(4-A) 9{Vii,-d,) = -p' 



Vp'(2-A) g{Wn.-d,). 



I.e., Vii and Wu are future pointing timelike in region II, and are candidate proper velocities 
of charged particles. Therefore, in region II we can write the current vector as the fiow of 
two oppositely charged, infalling timelike currents 
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where 



hi — hi a + hi b 



2,/p^{A-A) d_ r cosg 
a? sin 9 dO \_ sin^ 9 



Vn 



and 



2^p^{2 - A) _d r cosg 
sin 9 d9 \_ sin^ ^ 



C. Separating The Regions Outside The Event Horizon 

Collectively, h — I and /// = / give a meaningful description of the net current ev- 
erywhere in the magnetosphere. The regions we chose might lead us to believe that the 
difference in the character of our decomposition must occur at the event horizon. This is 
not necessarily the case. Vn and Wu continue to be future pointing timelike in some open 
interval outside the event horizon. I.e., there exists 5 > such that Vn and Wn are future 
pointing timelike in r < r+ + 5 (for a fixed the surface given by f = r+ is compact). 
Therefore in the region given by r < r+ + 5, we set 



I — hi — hi a + hib 



and for r > r+ + 5 we set 



I = h = ha + hb . 



In fact, the new regions has further advantages, in that 




la 



and 




lb 



are well defined. 



IV. EXTRACTION OF ENERGY FROM MATTER CURRENTS 



Here we focus on the outgoing current in region I: 



2^3/(3/ - 2) d r . cos 9 
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It is easy to see that Ijb is divergence free. We now use this conserved current, that is the 
only outflowing current in this decomposition, to construct an expression for the extraction 
of matter energy from the black hole. Define the charge density Pc by 

and the current 3- vector J by 

f = a Pj^ 

for i = 1,2,3. /jf, = imphes that 

dtPc + V-J = 0. 

Now, assuming that at every point lib is comprised of only one species of charged particle 
with charge q and mass m, we can define the mass density of the current as 

m 
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so that Pm > 0. This happens when 
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In the above equation, n is the outward pointing normal, meaning n = dr/ ^/g7r when 
r — !■ oo and n = —dr/^/g^ when r — )■ r+, and dA = ^J^eelw dip. Therefore, the matter 
extraction of energy £m from the black hole is given by 
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For our particular current Ijb, we get 
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which from eq.( !T6|) is greater than zero, thus allowing for extraction of matter energy from 
the black hole. Of course, the above equation is meaningful only when A is able to absorb 



the infinity produced by (sin^^) ^, which is an easy task since A is any arbitrary function 
of 9. To understand the colhmation effects, it is important to note that 
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Therefore, we see that (sin 6')"^ factor gives the currents a preferential polar jet like fea- 
ture. Here, unlike the case of the stationary electromagnetic fields, the mechanism by which 
charged particles carry energy from the black hole is apparent from the nature of the out- 
flowing currents. If indeed there are regions of antiparticle currents, we must make sure to 
modify the sign of q appropriately in that region. 



V. THE LORENTZ FACTOR AT LARGE DISTANCES 



Since the Boyer-Lindquist coordinates are asymptotically flat, from the t component of 
Wj in eg. (1151) . we see that the lorentz factor F of the ejected mass at infinity becomes 

3/ - 1 



^/3l{3l 



The above equation can be inverted to give 
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As r — !• 1 the energy colhmation factor / — oo, and when 1^2/3 the Lorentz factor 
r — i- oo. In particular, for any finite /, the ejected particles are so energetic that it never 
comes to rest even infinitely far away from the black hole. Clearly, the colhmation effects 
on r stem from our freedom in choosing a judicious I. The strength of the jets at the poles 
are also compensated by the intensity of the emitted particles. 



VI. SYMMETRY PROPERTIES OF THE FORCE-FREE EQUATIONS 

Consider a complete description of the fields and currents given by quantities p, J, E, D, 
B, and H. For each of these quantities will be define a dual object (the dual of a quantity 
A will be indicated by A) such that all the dual objects collectively describe a force-free, 
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stationary, axisymetric magnetosphere in Kerr geometry. The dual quantity will be very 
simply related to the original quantity, and yet, collectively, the physical content of the dual 
solutions will not be equivalent to the original one. The general features of this construction 
involve looking at the poloidal components of an object separately from the toroidal and 
the zeroeth component of the covariant formalism. This reason for this should be fairly 
clear: much like the background geometry, our assumptions require quantities to be time- 
independent (affecting the zeroeth component of a vector) and axisymetric (affecting the 
toroidal component of a vector). 

The dual charges and current are defined by 

p = p , Jt = Jt , Jp = -Jp. (18) 

Therefore if we define 

D = D, (19) 
we see that Gauss's theorem will be naturally satisfied since 

V-D = V- D = p = p. (20) 

For the other inhomogenous Maxwell's equation to hold, define 

Hp = Hp , Hx = — Hx . (21) 

Then 

(V X H)t = (V X Hp)t = (V X Hp)t = Jt = Jt, (22) 

and 

(V X H)p = (V X Ht)p = (V X -Ht)p = -Jp = Jp. (23) 

The second equality in the above equation holds only because the fields are axisymetric, for 
example 

(V X HY = e'-'^^ideH^ - d^He). (24) 
Therefore, for time-independent solutions it follows from eqs. ( l22l) and ( l23l) that 

-dtD + V xH = J. (25) 
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Now lets consider the homogenous Maxwell's equations. Having defined H and D, we have 
no more freedom in picking B and E. It is not difficult to see that 

Bp = Bp , Bt = — Bx , (2^) 

and 

E = E. (27) 

Clearly, the curl of E vanishes, and once again due to axis-symmetry, the divergence of 
B is trivial as well. Therefore, we have shown that the new quantities satisfy Maxwell's 
equations. It is just a matter of simple calculation to show that the dual fields and currents 
are force-free. Therefore, there exists a 3-vector u such that E = —u x B. It turns out that 
u = u. 

From eq.( l2T]) we see that if satisfies the Znajek regularity condition, will not 
(unless = —H^ = 0). Therefore, if we are using the dual solution to describe the external 
magnetosphere, we must separate the regions at r = r+ + 5 as explained in subsection IIII C[ 



VII. THE SOLUTION 



The dual solution to the f2_ solution presented earlier is given by 
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where / is the principle outgoing null geodesic of the Kerr geometry given by 

a 



+ a? 



We will refer to the set of equations above as the solution. 



A. Extraction Of Energy From The Electromagnetic Poynting Flux 



These solutions do allow for the extraction of energy via the electromagnetic Poynt 
flux from the black hole. From eq.f ll2p we see that the rate of energy extraction 



d£ 



EM 



dt 



j {-H^)nB'' ^rdA 



dS 



EM 



dt 



> 



For the case of our f2_ solution, the above expressions gives that 
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The collimation effects from the electromagnetic fields are given by 
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B. Extraction Of Energy From The Electromagnetic Currents 



In region I, the dual currents can be decomposed as follows: 
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Naturally, the extraction of matter energy stems from Ijb- Analogous to eq. ( !T7|) . here the 
matter extraction of energy Sm from the black hole is such that 
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which is positive when q{6) is correctly chosen so that 
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VIII. A PARTICULAR CHOICE OF LAMBDA 



Since we want energy extraction via both the electromagnetic fluxes and the matter 
currents, through out this section, we will focus on the Q- solution. Consider the simplest 
case where 

A = Ao sin^ e . (40) 

The factor of sin^ 6 is necessary to make the fields well defined on the poles. From eq. f l39p 
we have that 

3 71 

q = q- <0 when <9 < tt/A , and — < 6 < n . 

and 

q = q+>0 when 7r/4 < 6^ < 3 7r/4 

when Aq is positive, which we will now require. Let m_ and m_|_ be the corresponding mass 
of the particle species. Then 

d'^SM m 2 Ao cos 26 



dt dA 



q 0? sin 6' 



Clearly, this describes a jet like solution where the rate of energy extraction is maximized 
at the poles. There is a second local maxima along the equatorial plane suggesting the 
existence of a strong accretion disk. The above expression is integrated to give 
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We are not assuming that the species of charged particles in the two regions (near the poles 
and near the equatorial planes) are particle antiparticle pairs, although nothing precludes it 
in our formalism. They in could in particular be currents of electrons and protons (so that the 
black holes remains neutral during the process of energy extraction). The electromagnetic 
Poynting flux however gives that 

d'^^EM _ £ Aq cos^^ 

dAdt ^ ' ^ ' 

Here, just as in the case of the matter currents, the rate of electromagnetic energy extraction 
is a maximum along the polar axis. Interestingly, at the equatorial plane, the electromagnetic 
extraction rate is trivial suggesting a secondary mechanism for the observed glow around 
the accretion disk. The total extraction rate is given by 

d£EM 16 vr Aq 



dt 3 



(42) 



IX. MAGNETIC MONOPOLES 



Since W = , as long as A is the same in regions I and II, we will not introduce any 
magnetic monopoles as a result of using the dual solution in the external region I (r > r++ 5). 
However, it is important to pick A carefully if we are to exclude magnetic monopoles in every 
closed region of spacetime. Our choice of A given by eq. ( l40l) is one such choice. Indeed 
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X. CONCLUSION 

There are two important clarifications that we must make before we conclude this article. 
In we had claimed that it is impossible to extract energy when Q = Q^. The reference 
in this case is to the energy extracted via the electromagnetic flux alone; and it still holds 
true here (as it must). In this paper, as we have seen, it is possible to extract energy via 
matter currents when Q = Q^. The extraction of energy via electromagnetic flux also does 
occur when Q = Q_. The price to pay in this case is in the discontinuity of if^. This in turn 
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produces a delta function current at the membrane joining the two regions. The analysis of 
this current on the membrane requires further study. 

Outside of a few mild constraints, the functions I and A are arbitrary functions of ^. It 

is not clear whether astrophysical black holes permit a wide variety of magnetospheres, or 
if there is some other mechanism restricting the large degrees of freedom the decomposed 
currents and the fields have. Nonetheless, we have constructed a specific, exact solution to 
the Blanford-Znajek mechanism that extracts energy from the black hole. 
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Appendix A: Kerr Geometry Essentials 

For completeness, we define the various Kerr coordinates used. For asymptotic analysis, 
the Boyer-Lindquist coordinates are preferred, while the horizon and the interior region 
(r < r+) is analyzed using the usual Kerr-Schild coordinate system. 

1. Boyer-Lindquist Coordinates 

In the Boyer-Lindquist coordinates {t, r, 9, ip} of the Kerr geometry, the metric takes the 
form 



+7rr dr"^ + 'jee d6^ + 7<^<^ dip'^ , 



where the metric coefficients are given by 
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Here, 



and 



Additionally 



and 



p2 ^ + cos^ Q , A^r"^ - 2Mr + 
E2 = (r2 + a2)2-Aa2sin2^. 
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E2 



-\/^ = a = sin 9 . 

The parameters M and a are the mass and angular momentum per unit mass respectively 



of the Kerr black hole. The horizons are located at r± — M i: \J M"^ — a^. 



2. Kerr-Schild Coordinates 

Kerr-Schild coordinates are given by the transformation 
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In this frame, the metric becomes 
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where z = 2Mr/p^. Components of vectors transform as 



and 1-forms transform as 
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where the matrix A is given in eq.f lAip . 



We pick our time orientation for the Kerr geometry such that the null vector field —df is 
future pointing everywhere. 



Appendix B: Equations Of Electrodynamics In Stationary Spacetimes 



We only state the relevant equations of electrodynamics of stationary spacetimes. For a 
detail development, see 



Maxwell's equations can be written as 
* F"^ = , and V/jF"^ = 1° . 



(Bl) 



Here F"^ is the Maxwell stress tensor, J° is the four vector of the electric current and V is 
the covariant derivative of the geometry, -k F is the two form defined by 



(B2) 



Here, ea/Sfiu is the completely antisymmetric Levi-Civita tensor density of spacetime such 
that eoi23 = y/—g = <yVl (q^ ^ind 7 along with the other relevant Kerr quantities are defined 
in section |Al In the 3+1 formalism, where do is the asymptotically stationary timelike killing 
vector field, E and B are defined so that 
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We also define dual vectors D and H by 



*F, 






Hi 


H2 


Hs 


-Hi 









-H2 









-H3 










(B4) 



Naturally, F and * F are not independent. They are related by 

aD = E - (3 X B 



(B5) 
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and 



Here, 



H = aB - (3 X D . 



(B6) 



{A X By = e'^^ Aj Bk , (B7) 

where e^^'' is the Levi-civita tensor of our absolute space defined a;° = constant. Also, f3 
is the shift dual vector given by /3 = /3<^ d(p. Naturally, the spatial coordiantes are given 
by {x^,x'^,x^), and three vectors E, B, D, H live in this absolute space. Now, Maxwell's 
equations can be re-written as 

V ■ 5 = , (B8) 
dtB + V X E = , (B9) 
V-D = p,, (BIO) 



and 



-dtD + V X H = J 



(Bll) 



where pc = a/* and = al^. Here pc is the charge density and J is the electric 3-current. 
V is the covariant of the 3 space with the induced metric. The force-free condition that we 
will enforce is 

. (B12) 



^ VOL ^ 



This condition takes the form 



and 



J = 



PcE + J X B = 



(B13) 



(B14) 



For the case of a stationary, axis-symmetric, force-free magnetosphere, it is easy to show 
that there exists u = VL du^ such that 



E 



-uj X B . 



Additionally, [Uj showed that 



sin^ 
a 



(B15) 
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is the required condition in the Boyer-Lindquist coordinates that the otherwise bounded 
fields must satisfy so that they continue to be well defined in the Kerr-Schild coordinates at 
the event horizon. Eq. dBlSll is referred to as the Znajek regularity condition. 
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